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CHARGE-EQUILIBRIUM AND RADIATION OF LOW-ENERGY 


COSMIC RAYS PASSING THROUGH INTERSTELLAR MEDIUM 

D, W. Rule* and K, Omidvar 
Stratosphere Physics and Chemistry Branch 
Laboratory for Planetary Atmospheres 
NASA/Goddard Space Flight Center 

ABSTRACT 

The charge equilibrium and radiation of an oxygen and an iron beam in the 
MeV per nucleon energy range, rep resenting, a typical beam of low-energy cos- 
mic rays passing through the interstellar medium, is considered. Electron loss 
of the beam has been taken into account by means of the first Born approximation 
allowing for the target atom to remain unexcited, or to be excited to all possible 

r 

states. Electron capture cross sections have been calculated by means of the 
scaled Oppenheimer-Brinkman-Kramers approximation, taking into account all 
atomic shells of the target atoms. Radiation of the beam due to electron capture 
into the excited states of the ion, collisional excitation and collisional inner- 
shell ionization of the ions has been considered. Effective X-ray production 
cross sections and multiplicities for the most energetic X-ray lines emitted by 
the Fe and O beams have been calculated. 

Subject headings : atonic and molecular processes - cosmic rays - 
interstellar matter - X-rays. 

*National Research Council Research Associate 
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chaege-equilibriiim; and radiation of low-energy 
COSMIC rays passing through interstellae medium 

I. INTRODUCTEON 

When a neutral or charged beam of particles passes through a gaseous 
medium, through collisions, it reaches a state of charge equilibrium independent 
of the initial charge state of the beam. The equilibrium is reached through the 
competing processes of ionization of the beam by the surrounding gas, and elec- 
tron capture by the beam from the gas, the equilibrium time being of the order 
of the collision time. 

Through the colllsional excitation of its particles with the surrounding gas 
the beam radiates. The radiation is in die visible to the X-ray range, depending 
on the composition of the beam, is Doppler shifted in frequency if the beam is 
unidirectional, and is Doppler broadened if the beam has an isotropic distribution. 

In a number of problems of interest such as the propagation of isotropic low 
energy cosmic rays through the interstellar medium, or neutral beam injection 
in a fusion plasma, information about the charge state equilibrium of a passing 
beam and its radiation distribution and intensity is desirable. This information 
requires an accurate knowledge of cross sections for charge-exchange, ioni- 
zation, and excitation processes. 

There have been several calculations of charge equilibrium and X-ray pro- 
duction by oic^gen and Iron ions, e.g. , Serlemitsos 1973; Watson 1975; Pravdo 
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and Boldt 1975; and Bussard, et al. 1978, However the previous work contains 
a variety of approximations for the charge exchange, ionization, and excitation 
cross sections needed as input, and consequently the results of these calculations 
are not in agreement in several instances. Recently the reliability of the iron 
ion calculations using scaled Brinkman-Kramers charge exchange cross sections 
has been questioned (Belki8 and McCarroll 1977). For these reasons, the 
present work focuses on the method of obtaining the input cross sections, all of 
which have been consistently calculated within the Born approximation. By 
careful comparison with experimental data wherever possible, we have attempted 
to assess the reliability of the methods we have used. In this way we have been 
able to understand the sources of some of the disagreement among previous work. 
Furthermore, our analysis should provide a reliable' starting point from which 
to approach similar calculations involving different atomic species. 

In the next section a brief description of the theory will be given. In the 
third and fourth sections results and discussions will be presented. 
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n. CALCULATEON OF CROSS SECTIONS 
a) Ionization Cross SectioKS 

The ionization cross sections used in our determination of the equilibrium 
distribution of charge states were obtained using the first Born approximation 
(FBA). These calculations are an outgrowth of the methods developed for hy- 
drogen and helium projectiles by Bates and several co-workers (for a review 
see Bates 1962). Similar calculations of electron loss by hydrogen-like and 
helium-like ions in neutral hydrogen and helium have been done by Dmitriev, 
et al. (1965) and Senashenko, et al. (1968). Recently, this approach was ex- 
tended to the calculation of ionization from multi-electron atoms in collision 
with neutral atoms by Rule (1977), In this section we will simply outline the 

calculational method. The references given in this paragraph contain more 

1 

detailed discussions. 

We begin by assuming that the major contribution to electron loss by an ion 
is direct Coulomb ionization of a single electron. We thus neglect Auger pro- 
cesses, multiple electron ionization, and transfer of electrons from the pro- 
jectile to any target ions that may be present. We then write the total ionization 
cross section as a sum over partial cross sections for each occupied subshell 
of the projectile ion: 

N 

»'(V) = (1) 

n = 1 2 
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where N denotes the principle quantum number of the highest occupied shell, 
and V is the relative velocity of the ion-atom system. Each g can be sep- 
arated into a part which represents ionization processes in which the target is 
elastically scattered, aj. , and a part which represents the sum over all ioni- 
zations of the projectile ion during which the target' atom is inelastically 
scattered, o\ , as 

^n,£W - + o\m (2) 

The FBA expression for these cross sections is (see, e.g. , Bates 1962): 

Zt- /*^2 dQ 

ui(V) = 87ra2 -p— / _ li-F(Zn^gQ)P / k?dk8.^^^(k,Q), (3) 

°(Z2gv)2-/Q, Q3 Jo 

and 

z r'^2 do /■''* 

o[(V) = 8ira2 _ S(Z„,(.Q) k2d)c&„,«(k,Q). (4) 

In the above expressions, slq is the Bohr radius, g is the effective charge for 
an electron in the (n, £) subshell of the projectile, and 

V s V/Z„,JV„ (5) 

where Vg is the Bohr velocity. We also have Z-p, the target's nuclear charge, 
and Q, which is equal to times the momentum transfer, 

Q = Zn‘,E(Kf-Ki>. (8) 

— y 

Here Kj and Kf are the initial and final momenta in the center of mass of the 
projectile-target system in units of a^^ . 


4 



7 


The functions in the integrands above are F(q) , the target elastic' form 
factor for momentum transfer and S(q), the Incoherent scattering form factor 
for the target atom. The form factor for transition from the (n, 2) subshell to 
the continuum state with momentum of magnitude k is given by 

N 0 ^ 

^n,eCk,Q) = dk Y, (2-5^,o)l<klei^*?ln2m>|2, (7) 

22 + 1 y = 0 

where q is the number of electrons in the subshell and I k) and In 2m) are 
hydrogenic continuum and bound state wave functions. The explicit form for 
^n,C (k»Q) Is given by Omidvar, et al. (1972). 


The limits in Equations (4) and (5) are as follows: 

Q, = In.c/aZ^S'') • 
Q„Qi « 

and 

Q; CIn,2 + Ix)/(2Z2^gv). 


( 8 ) 

(9) 

( 10 ) 


The upper limits of the k integrations are given by 


ki = 


mp 


L Pt 


-1 1/2 


(2QKi - Q^) - In,2 


( 11 ) 


and 


k! 




1 


n 1/2 


(2QK; - Q2 )--(I „,8 + It) 
L 2 m 2 


( 12 ) 


where mg is the mass of an electron, and ju is the reduced mass of projectile- 
target system. Also, is the ionization potential for an electron in the (n, 2) 
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sub shell in Rydbergs and is the ionization potential for the neutral target 

atom. Note that when v> 1, the proportionalities 

(7e(v) CC and oc (13) 

give scaling with respect to Z^ji and v, which is identical for both crj and aj, 
while the scaling with differs in the two cases. 

Ionization cross sections for Pe'^20 hydrogen are shown in Figure 1, 

The solid curves were calculated as described above, while the dashed curves 
were obtained from the binary encounter approximation (BE A) as tabulated by 
Hanson (1973) for ionization by protons. The data of Berkner, et al. (1977) is 
shown for comparison. Since the original data was taken for F e-ions in IJ 2 , we 
have divided it by two before plotting. We see that the agreement of the present 
calculation with the data is quite good. For the case of Fe"^^^ and Fe"^^^ , the data 
falls between the BEA and the FBA values; however only a single data point is 
available for these cases. 

We have found that at energies near the peak of the cross sections, the 
contributions of OgUnd aj to Equation (2), are comparable. For example, 

at ~25 MeV/nucleon R = a|/ag==^0,6. At energies lower than that of the peak, E 
becomes smaller, e.g. , at 0, 67 MeV/nucleon R 0. 5. The contribution of o\. 
is roughly equivalent to the cross section for ionization of Fe-ions by a proton 
(see Rule 1977) as was used by Watson (1975), thus neglect of the structure of 
the hydrogen target atom during the ionizing collision can cause a substantial 
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error in the calculated cross sections. ]For lower charge states of T'e, - 

the ratio ~0*5 near the peak of the cross sections. The ioni- 

zation cross sections for the ions Fe"^^® - Fe’*'^® have also been calculated in 
the same manner as those in Figure Fusing the FJBA. ■ - 

In Figure 2 the calculated ionization cross sections for single electron loss 
by oxygen ions are given. The FBA results are again compared to the BE A 
results. We are not aware of any data for comparison with these calculated 
values; however, for the ionization of O'^'^ - Q'*'^ by He, a detailed comparison of 
the FBA and BEA values with the data of Macdonald and Martin (1971) has been 
given by Rule (1977), There it was shown that FBA results were in reasonable 
agreement with the data for all five charge states considered, while the BEA 
values over-estimated the cross section considerably for the O'*'^ - O’*'^ cases and 
underestimated them for and The relative contributions to £ from 

Og and a I of Equation (2) are again indicated by the ratio E = " 

O’*'^ in H, R « 0. 5 near the peaks of these cross sections, while for in 

He E 0. 3 for peak values. Thus inelastic target processes collectively con- 
tribute substantially to the total ionization cross sections for oxygen ions as well 
as for iron ions. 

b) Excitation Cross Sections 

The first Born approximation approach was also used to calculate the ex- 
citation cross sections necessary to obtain effective X-ray production cross 
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sections. The method is quite similar to that outlined above for the case of 
ionization. The ealculational technique has been described thoroughly by Bates 
(1962) and also by Levy (1969). We vidll simply give the results here. 


As for Ionization, we can divide the total excitation cross section into 
two partsj .Og and a?, which represent excitation of the projectile during elastic 
and inelastic scattering of the target, respectively. Just as in Equations (3) and 
(4) we can write: 





dQ 

F(Z„,jQ)|2 — , 


and 

Zt dQ 

The function ^nK,n'£'(Q) defined in a manner analogous to g(k,Q) of 
Equation (7): 


(14) 


(15) 


^n'2,n'£'(Q) ~ ^m,o) ^ ^ 


( 22 + 1 ) 


m = 0 


As in Equation (7), InSm) is a hydrogenic wave function. 


gration in Equations (14) and (15) are: 


ei^’^ln2m)P, (16) 
The limits of inte- 


and 


Q2.Q2 « 

Qi = Qi + L,/(2Z2^gv), 


(17) 

(18) 


(19) 


8 



11 


where is an effective excitation energy for the target atom. For excitation 
of light atoms such as H, Q'^ is sensitive to the value of I-r (Levy, 1969), which 
may be taken as the average excitation potential determined in the Bethe-theory 
of stopping power (see Dalgamo (1962) and Fano (1963) for a review). For H, 

Ij = 15, OeV; for He, = 41.8eV, For highly charged heavy ions, It is appar- 
ent that the first term in Equation (19) will dominate the second; therefore, the 
expression for Q'j^ will not be sensitive to the value of Ij. This was found to be 
especially true for the excitation of Fe-ions. 


For highly charged ions, ^ 1, then > I^/(2Zn^g v) and «=* Q^, also 
if V > 1 Q'^ 0. In either case, the approximate equality of the limits of 

integration in Equations (14) and (15) allows the integrands to be combined to 
give 


(jE(v) = (jE + g^a2 



dQ 

— [1 -2F(Zn {Q)] • S-n£,n'S'(Q)‘ (20) 

Q 


For hydrogen targets Equation (14) becomes (Bates, 1962) 

1 /*“ dQ 

ffE(v) ft:* 167ra2 — / — [1 -F(Zn,£Q)] Snfi,n'£'(Q)> (21) 

which is roughly a factor of two larger than o|(v) of Equation (15), For v 1, 
and for > 1, F(Zji^g Q^) ^ 0, i. e. , for highly charged ions, the minimum 
momentum transfer for ionization is so large that, at that momentum, the 
target form factor is small and does not contribute significantly to the integral 
in Equation (21), 
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For Fe-ions in atomic hydrogen we have found Equation (21) to be quite 
adequate for calculating excitation cross sections; however, in the case of 0-ions, 
we found that the approximations made in obtaining Equation (21) were not jus- 
tified, and we calculated cr| and up separately. Here we should point out that 
excitation by protons is approximately given by u|(v) alone. In the case of ex- 
citation of O-ions, the present results are roughly a factor of two larger than 
those of Pravdo and Boldt (1975) which were calculated using only the formula 
for uf(v). Equation (14), 

Cascading from higher excited states to the upper level of the observed 
transition was also Included in our calculation. For example, the total (Is 2p) 
excitation cross section was calculated as 

nc- n - 1 

|3T(2p,nS) [af(ls-^n,£) + af (ls-^n,£)] , (22) 

n = 2 C = 0 

where /3Y(2p,nJ2) is the total branching ratio for (n2 2p) transitions (Omidvar 

1978), and where Uc = 6, contributions for n > 6 being negligible. 

Figure 3 contains examples of cttotal(^®~^ ^P) calculated as described 
above. Cascading contributed ~20% to the - H and - He cross sections 
at 0. 1 MeV/nucleon, while at l.OMeV/nucleon cascading contributed ~10% in 
each case. For Is ^ 2p transition of Fe"^^^ in H (not shown), cascading con- 
tributed ~20% to Ototal(^® ^ ^P) 1 MeV/nucleon and ~10% at 10 MeV/nucleon, 

which is at the peak of the cross section. 
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c) Charge -Exchange Cross Sections 

In the passage of charged particles through gases, the main process for 
neutralizing the charged particles is the capture of electrons from the surrounding 
gas. If the surrounding gas is partially ionized, two other processes, dielec- 
tronic recombination and radiative recombination, can also contribute to the 
neutralization of the beam. The dielectronic recombination which is the inverse 
of the Auger effect can be neglected if the beam is made up basically of bare 
nuclei. Radiative recombination which is equivalent to the free-bound spon- 
taneous emission has a cross section which is approximately a factor of 
smaller than the charge-exchange cross section, where a is the fine structure 
constant, 

A fourth process which can occur with neutral surrounding gas is the radia- 
tive charge-exchange process (Briggs, et al. 1977; Butler, et al. 1977) but 
this is unimportant at energies of our interest. 

There are numerous calculational methods applied to charge-exchange 
collisions, but they are mostly limited to charge exchange between protons and 
the simplest atoms, such as atomic hydrogen or helium. For highly charged 
ions charge-exchanging with multielectron atoms in the MeV energy range, 
where other than the K-shell electrons contribute substantially to the charge- 
exchange cross sections, calculational methods are scarce. Aside from the 
Oppenheimer-Brinkman-Kramers approximation that will be discussed below. 
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two other approximations can readily be applied to the multielectron case. One 
is the method of Jackson and Schlff for electron capture by protons from atomic 
hydrogen, where the interaction between the proton and the H nucleus is taken 
into account (Jackson and Schiff 1953). This method, with some modifications 
when applied to the K-electron capture from He and Ar, give satisfactory agree- 
ment with measurements (Omidvar, et al. 1976). 

The other method that is being considered for application is the Bates' 
Atomic Expansion Method (Bates 1958), An Important feature of this method 
is that it accounts for the lack of orthogonality of the initial and final atomic 
states. This method has been applied by Lin, et al, (1977) to the argon K- 
electron capture by protons. Good agreement with measurement for both the 
total and the differential cross section is obtained. 

The calculational method used here is the Oppenheimer-Brinkman-Kramers 
approximation (Mapleton 1972) modified by an empirical factor (Nikolaev 1967), 

Let a{xi 2^2 > “■ 1 ^ 1 ) represent the charge exchange cross section for capture 
of an electron that is initially in a state of principal and azimuthal quantum num- 
ber vi 2^2 of fh® target atom, and finally in a state n^ f i of the projectile. The 
cross section for population of the state Uq of the projectile is then given by 

Ox ( 02 ^ 2 .^ 020 )= ^ (23) 
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where |S-p(ni£i,noJZo) are the electric dipole branching ratios for nQ^o 

transitions. The a(n 2 fi 2 jaiSi)'were calculated using the analytic form of the 
charge exchange amplitude for the hydrogenic wave functions given elsewhere 
(Omidvar 1975). Extensive tabulation of is^ given by Omidvar (1978). 

j 

In Table 1 some recently measured cross sections are compared with related 
calculated cross sections. As is seen, the calculated charge- exchange cross 
sections are in fair agreement with measurement for H and He targets, but for 
atoms that have other than K-shell electrons, the calculated cross sections are 
too large by a factor of 2-6, and in some exceptional cases this factor is as 
large as 23. In Section IV we will make an estimate of the change in our results 
if we reduce the charge exchange cross section by a factor of 5. 

The charge exchange cross sections are presented in Figures 4 and 5. 

In a number of ways the present charge-exchange calculation is an improve- 
ment over the previous calculations (Seriemitsos, et al. 1973; Pravdo, et al, 
1975; and Watson 1975), The most important improvement is to account for the 
contribution to the charge exchange cross section by the target atomic shells 
other than the K-shell, while these authors have consistently negelected the 
non -K-shell electrons. From Figure 4 it is seen that the non-K-shell contri- 
bution to the cross section is substantial and in some cases dominate. An 
account of other improvements can be found in Bussard, et al. (1978). 
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III. CHARGE EQUILIBRIUM AND -RADIATION 


a) Charge Eractions 


The ionization and charge exchange cross sections given in the last sec^n 
were used to calculate the fraction of the ion-beam in each charge state as a 
function of energy,. The rate of change of the fraction F 2 of. a beam with charge 
Z is given by (e.g. , see Nikolaev 1965) 

dF, 


= E Pj <'jz - Pz Z "Zj , 

dn j z j z 


(24) 


where du = pdx, with units of atoms/cm^. Here dx is an interval of distance 
over which the velocity of the projectile can be considered constant and p is the 
density of target atoms in atoms/cm^ . The cross section Ojz is the probability 
of an ion initially of charge j having a final charge Z after interacting with a 
target atom. As dF/du = 0, and an equilibrixim charge distribution is 

established in the beam at each velocity. Under these conditions we can express 
the condition of equilibrium as 

Z Pj = Z Fj > (25) 

J < z j > z 

K > Z K < Z 

If we assume an equilibrium charge distribution exists for low energy cosmic 
rays, then we have an equation like Equation (25) for each charge state Z, By 
further assuming that charge changing processes involving only one electron are 
dominant, Equation (25) reduces to 

Fz ~ F2+1 <^z+i > ( 2 ®) 
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with ^ F 2 = 1 and where 
z 

= ^z,z + i “ ‘^z+i.z (27) 

are the single electron ionization and charge exchange cross sections* 

In Figure 6 we show the result of determining the charge fractions from 
Equation (26) for a beam of 0-ions in atomic hydrogen for charge states to 

Also shown is a similar result obtained for a mixture of gases representing 
the interstellar mediuna. For this case, effective ionization cross sections per 
H atom were calculated as 

+ 0.085 a^CHe), (28) 

where the helium cross section is multiplied by its abundance in the interstellar 
medium relative to hydrogen. The effective charge exchange cross sections per 
H atom were taken to be 

A Nh 

i.e. , the sum over the values for the most abundant elements times their abun- 
dances relative to hydrogen. The elements contained in the sum are those in- 
dicated in Figure 5. Since cr^ scales roughly as the square of the target atom's 
nuclear charge, , while 0 ^ scales roughly as zf- to Zj , less abundant ele- 
ments contribute significantly to , as discussed in Section II, while only He 
contributes significantly to o^, Equation (28). Also shown for comparison are 
the set of charge fractions determined by Pravdo and Boldt (1975). The effect 
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of their use of charge exchange cross sections which are too large by a factor 
of about 3 is to slightly shift the charge fractions toward higher energy. 

In- Figure- -7- the- set of charge fractions- for a- beam-of -Fe-ions- in -H atoms 
with charge states from Fe"^^^ to Fe^^® are given. Also shown are the charge 
fractions for the case of the interstellar mixture of gases, M this case 
was determined by scaling the elastie-and inelastic-target cross sections [Equa- 
tions (2)-(4)] separately as; 

"z = S — 4®) + Za oJCH)] . (SO) 

Na 

while the charge exchange cross sections were obtained as in Equations (23) and 
(29), The effect of the interstellar medium was to shift the charge fractions to 
higher energies by ~3MeV/amu, and also to broaden the charge fraction curves. 
This broadening seems to be caused by the steeper slope of in pure H atoms 
as compared to in the interstellar mixture (see Fig, 5), 

b) Effective X-Ray Production Cross Sections 
The effective X-ray production cross sections were calculated by multi- 
plying the appropriate cross sections for excitation, ionization and charge ex- 
change by the fractions F 2 of the beam having a given charge state. For ex- 
ample, to calculate ax(0.65), the cross section for ojgrgen ions to produce 
0. 65 KeV X-rays, we take 

Ox(0.65) = Fg 08^(2p) + F, of (1 s ^ 2p), (31) 
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where a^( 2 p) is for capture of an electron by into the 2 p state, including 
cascading from capture into higher levels. 

In Figure 8 we have plotted FgOg (2p), and (ls^2p) for 

X-ray production by O-ions in the interstellar medium. There it is seen that 
the contributions from charge exchange and excitation are nearly equal. This 
result seems to be in accord with that of Watson (1975), while it disagrees with 
the Pravdo and Boldt (1975) result that charge exchange dominates excitation by 
roughly an order of magnitude. The source of this discrepancy is probably their 
underestimation of a® (is - 2p) by a factor of two, as described In Section II , and 
the overestimation of Ug (2p) caused by using unsealed Brihkman-Kramers cross 
sections, which are a factor of about 4 too large. 

Figure 9 is a plot of the cross section for production of ~0.57KeV X-rays by 
O-ions in the interstellar medium. In this case we have taken 

ax(~0.57) = F 7 ff^(2p) + Fg af (ls-2p) (32) 

Here we have neglected the contribution to f^x(~0,57) from inner shell ionization 
since the K-shell fluorescence yield Wjj; is small. For example, for neutral 
o:?iygen ** 0,009 (Bambynek, et al, 1972), We point out that the expression 
for in Equation (32) differs from that of Pravdo and Boldt (1975) by the in- 
clusion of the last term, while it is identical to that of Watson (1975). 

In Figure 10 we give the results for X-ray production by O-ions in atomic 
hydrogen for comparison with those of Figures 8 and 9 obtained with the 
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interstellar medium. In Figure 10 both ajj(0.65) and 0. 57) are plotted. 

The contributions to ax(~0.57) are as described by Equation (32). 

-Figure 11 -is. a-plot.o£ %-(6..9-7.),.-the..effectLv-e-Xrray_produc.tion cross, section, 
for Fe-ions in atomic hydrogen (labeled H) and in the interstellar medium (labeled 
I). In both cases is given by 

ffx(6.97) = F 26 u^(2p) + F 25 of5(ls^2p), (33) 

where the notation is the same as in Equation (31). Although not shown, the 
contributions of the two terms in Equation (33) are roughly equal, as in the case 
of ojq'^gen shown in Figure 8 , 

Figure 12 is a plot of Ux(~ 6 , 68 ), for X-ray production by Fe-ions in the 
interstellar medium, 

a^(~ 6 . 68 ) « F 25 a^5(2p) + F 24 a;^(ls-^2p) + co^ S ^z- (34) 

z > 22 

In Equation (34) we have included a term which gives an estimate of the effect 
of inner shell ionization on ctx(~ 6 . 68 ). This estimate could be improved by re- 
placing (Is^) by the true inner shell ionization cross section for each charge 
state, and by replacing the fluorescence yield for neutral Fe, , by the values 
corresponding to each charge state co| . Thus the last term in Equation (35) 
should be 

y ' Fz (35) 

z > 22 
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however the quantities are, to our knowledge, not available in the literature. 
For this reason our estimate in Equation (34) uses 0.437 which is an aver- 
age yield for Fe'*'^^ - Fe'*'^^. A more complete discussion of inner shell ioniza- 
tion contributions is given in the Appendix. 


Figure 13 gives ~ 6 . 68 ) for X-ray production by Pe-ions in atomic 
hydrogen for comparison with Figure 12. 


In Figure 14a and b we have plotted the effective X-ray production cross 
section for the 3p - 5 - 2s transitions of Fe at ~1. 13 to 1. 17KeV and also for the 
3d 2p and 3s 2p transitions at ~ 0 . 7 to l.lKeV. Figure 14a was obtained 
for the interstellar medium and Figure 14b for atomic hydrogen. These cross 
sections were obtained from the expressions 

(7x(3p^2s) F 23 a|3(2s^3p) + F 24 [ 0 ^ 4 ( 25 ) + 0 ^ 4 ( 2 p)] 

(36) 


+ F 22 022(2s“^3p) "I" F 23 [c7^(2s) + o^( 2 p)] , 


and 


21 


22 


ax(3d,3s-^2p) X] F^ af(2p-»3s,3d) + Fz a§(2p). 


(37) 


16 


Z = 17 


As in the previous cases, the interstellar mixture has the effect of broadening 
the energy range in which is significant. Here we have neglected contributions 
from inner shell ionization of the 2 p level, which should be small. 


We should point out that, in Equations (36) and (37), the cross sections for 
capture by Pe+^i - Fe"^!^ were obtained by scaling with respect to Zjj g as, e. g. , 
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( Z \ ^ 

(7^3 (2 p), (38) 

^ 2 , 0 / 

where the a's for each case were determined from the calculated cross sections 
for a^'(2p)' and' o^(2p) for each energy. For this reason Ox (3p^'2's)‘ and <7x(3d, 
3s -^2p) should not be as reliable as the other effective X-ray production cross 
sections. 


c) Multiplicities 

Table 2 is a summary of results for the calculation of the total number of 
X-rays produced by a cosmic ray nucleus as it slows down in atomic hydrogen 
or in the interstellar medium. For a given X-ray production process, the 
multiplicity, M, is calculated as 

^ dE, (39) 

where and Ef are the initial and final ion energies in MeV/amu chosen such 
that ffx (E) ^ 0 outside the range Ef < E < Ej . The effective X-ray production 
cross sections Ux cm^/(H~atom) have been described in detail above} and 
-dE/dX Is the stopping power of the medium in cm^-(MeV/amu)/(H-atom), which 
is taken from Northcliffe and Schilling (1970) and from Barkas and Berger (1964). 

Table 2 also contains results of calculations by Watson (1975), Pravdo and 
Eoldt (1975), and Bussard, et al, (1977), In comparing the present results with 
those of previous workers, it should be noted that the processes contributing to 
Ox differ among the various results cited. For example, in comparing the 
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present calculation to that of Watson for an atomic hydrogen medium, we find 
that, in each case, the latter multiplicities exceed our values. This is probably 
a result of Watson’s larger charge exchange cross sections and contradicts the 
conclusion of Belkie and McCarroll that Watson's results may be seriously under - 
estimating the X-ray production rate by and Fe'*‘25, in comparing our O-ion 

multiplicities with those of Pravdo and Boldt (1975), we find reasonable agree- 
ment; however this is misleading since the contribution to cr^ from charge ex- 
change and excitation is quite different. The Pravdo-Boldt was dominated 
by the charge exchange contribution, which exceeded excitation by an order of 
magnitude (see Section II), while we find the two contributions to be comparable, 
in agreement with Watson. Our results for Fe-ions in the interstellar medium 
are in reasonable agreement with those of Bussard, et al. (1977), with respect 
to both the multiplicities and the nature of the cross section Ox , whose largest 
component comes from excitation of the ion, rather than charge exchange. 
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IV. DISCUSSION 

In this section we will consider the effect which overestimating charge ex- 
change cross sections has on our final results for ion beams in the interstellar 
medium. In addition, we discuss the effect of multieleetron charge changing 
processes on the calculation of the equilibrium charge distribution. Finally, 
the contribution to X-ray production from inner-shell ionization of neutral Fe 
in the interstellar medium by low energy cosmic ray ions will be discussed. 

For the case of Fe-ions in atomic hydrogen, we have been able to compare 
our results for both ionization and charge exchange with the data of Berkner, 
et al. (1977). Agreement between calculated and experimental values was 
found to be within ~25% for hydrogen targets. Similar agreement was also 
found for oxygen ions in helium. However, comparison of charge exchange 
cross sections for Fe-ions in argon indicates that the Brinkman-Kramers approx- 
imation scaled with the Nikolaev factor does not correctly reproduce the empir- 
ical values. Thus our effective X-ray production cross sections for Fe-ions in 
the interstellar medium are not as reliable as they are for Fe-ions in hydrogen, 
since, for the former case, capture from less abundant atomic Fe and O was 
found to dominate. 

We have estimated the effect of using charge exchange cross sections which 
may be a factor of -^5 too large for Fe-ions in the interstellar medium. In 
Figure 15 we give the results of calculating ax(~6. 68) (Eq. 34) after dividing the 
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charge exchange cross sections by 5, Thus the equilibrium charge fractions 
and capture cross sections of Equation (34) will be altered from the values used 
to obtain Figure 12. The multiplicity for Ox(~6,68) from Figure 15 is 5.2 
which compares with 9. 8 from Figure 12. Similar effects should occur for 
the case of oxygen ions in the interstellar mixture. 

Aside from possible errors in the values of the single electron ionization 
and charge exchange cross sections which were used to calculate the equilibrium 
charge fractions and X-ray production cross sections, the processes involving 
multiple capture and loss of electrons may be important for the case of heavy 
targets. The data of Macdonald and Martin (1965) show that for O'*'® in argon at 
2. 5MeV/amu, double capture is about 15-20% of single capture. For Fe'''^® - 
-pQ +26 jn argon, in the range 1. 1-8. 4MeV/amu, double capture is ~20-30% of 
single capture, according to the data of Berkner, et al. (1977). For 
in argon at 1,1-2, 5MeV/amu, double loss is ~10 ~25% of single loss. For 
lighter targets these effects are much smaller. 

If one includes both one- and two-electron charge changing processes in 
determining the equilibrium charge distribution, then Equation (25) reduces to 
the set of equations 

^z-i^z-i,z+i ■^Pz[<^z,z+i ■‘■^z.z+zl ~ ^z+it^z+i.z <^z+i ,z-i 3 +Fz+ 2 °^z+ 2 ,z , (^0) 
one for each charge state Z. These equations can be put into a form similar to 
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Equation (26), except that and o^+i each become functions of several capture 
and loss cross sections: 

= ^ i ■ (41) 

z 

We have found that the equilibrium charge fractions based on single electron loss 
are relatively insensitive to changes in the values of (Eq. 26), as 

can be seen from Figures 6 and 7. Thus the reduction of Equations (40) to the 
form (41) suggests that as long as two electron charge changing processes are 
of the order of the one electron ones, the effect on the charge distribution of the 
beam should be small. However, we point out that when the target nuclear 
charge is greater than that of the projectile, multiple electron processes may 
dominate the one- electron ones and the above argument will then break down. 
Calculations of cross sections for multiple electron processes have been con- 
sidered, for example, by Bottcher (1977), 

Our crude estimate of the contribution to X-ray production by inner shell 
ionization of Fe-ions indicates that this process does not significantly increase 
the effective production cross section; however a consideration of Inner-shell . 
processes for the Fe-ions leads one to examine inner shell Ionization of the 
neutral Fe in the Interstellar medium. Although the abundance of Fe is ~10"5 
relative to hydrogen, the flux of protons and heavier ions which can ionize the 
K and L shells of neutral Fe is much greater than the flux of Fe-ions. The 
doppler broadening of this line will vary with the mass of the projectile. For 
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ionization of Fe by protons, the line should be quite narrow, as opposed to the 
result of heavy particles in, head on collisions. 

The K-shell vacancy production rate is proportional to the sum of ionization 
and capture cross sections (Halpern and Law 1973, LapicM and Losonsl?y 1977): 

a = cr^ + (42) 

For protons, is of order at its peak value, while is of order 

10'^ ^cm^ (see the curve labeled Is in Fig. 1), Thus a ^ for protons; how- 
ever rapidly increases for larger projectile charges, eventually dominating 
, e.g, , for O’*"® in argon, (Halpem and Law 1973). The values of 

can be obtained for heavier ions by scaling the results given by Rule for ioni- 
zation of Fe+25 protons. The appropriate a^’s can be determined from em- 
pirically scaled Brinkman-Kramers cross sections as discussed in Section II, 
Thus all of the data necessary to calculate the inner-shell ionization cross 
sections for Pe atoms by low energy cosmic ray ions are readily available. The 
effective X-ray production cross section is then 

Ox = (43) 

where is the K-shell fluorescence yield. A similar expression holds for 
X-rays resulting from L-shell vacancy production. 
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V. CONCLUSION 

We have used improved cross sections for charge exchange, ionization, and 
excitation of ions by neutral atoms, to calculate the charge equilibrium and radi- 
ation of an atomic beam passing through both atomic hydrogen and the inter- 
stellar medium. The present calculation results in a shift -with respect to energy 
in the distribution of the charge fractions of an oxygen beam passing through 
Interstellar matter. The magnitudes of different charge fractions are similarly 
changed. 

We have also calculated the munber of X-rays emitted per cosmic ray 
nucleus. We find charge exchange and excitation contribute comparably to X-ray 
production, while contributions from inner shell ionization are small. Partial 
agreement was obtained with the results of some authors; however, in some 
cases discrepancies were found. For example, we find our results for X-ray 
production by Fe"^24 ^nd in hydrogen to be about 1/3 the value obtained 

previously by Watson, 

Further Improvement in the charge equilibrium and radiation calculation 
necessitates a more accurate knowledge of charge -exchange cross sections, 
especially for cases involving capture of L- and M- shell electrons, 

^ 4 : 
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Table 1 

Comparison of some experimental and calculated cross sections and their ratios for a num- 
ber of charge-exchange reactions. E and are the energies and charge numbers of the 
projectile. o(Ex) and a(CAL) are the experimental and calculated cross sections in units 
of cm^. The first three cases of charge exchange are due to Berkner, et al. (1977), and 
the last case is due to Macdonald and Martin (1971). Random experimental uncertainties 
for the first three cases are approximately ±5%, and systematic uncertainties for these 
cases are ± 1 %. For the last case the experimental error is ~10-15%. 




+ H ^ 

+ H* 


E(MeV/amu) 


ct(Ex) 

a(CAL) 

a(CAL)/a(Ex) 

3.4 

25 

3.9 -20* 

3.49-20 

0.89 

3.4 

24 

3.6 -20 

2.71-20 

0.75 

3.4 

23 

3.0 -20 

2.36-20 

0.79 

3.4 

22 

2.7 -20 

1.98-20 

0.73 

3.4 

21 

2.3 -20 

1.71-20 

0.74 

3.4 

20 

2.2 -20 

1.47-20 

0.67 

1.1 

21 

5.0 -18 

4.72-18 

0.94 

1.1 

20 

4.0 -18 

4.13-18 

1.0 



+ At ^ Fe'^^^i 

+ Ar" 


E(MeV/amu) 

Zi 

(j(Ex) 

or(CAL) 

a(CAL)/a(Ex) 

8.4 

26 

5.8 -18 

1.62-17 

2.8 

8.4 

25 

5.6 -18 

1.44-17 

2.6 

3.4 

25 

4.5 -17 

2.57-16 

5.7 

3.4 

24 

4.1 -17 

1.91-16 

4.7 

3.4 

23 

3.7 -17 

1.67-16 

4.5 

1.1 

21 

1.46-16 

3.42-15 

23.0 

1.1 

20 

1.43-16 

2.94-15 

21.0 


Ar-^^i 

+ At ^ Ar"^^^! 

+ Ar^ 


E(MeV/amu) 

Zi 

a (Ex) 

or(CAL) 

a(CAL)/(r(Ex) 

8.4 

18 

2.4 -18 

5.78-18 

2.4 

8.4 

17 

2.4 -18 

4.96-18 

2.1 

3.4 

16 

2.3 -17 

5.35-17 

2.3 


Q-’-Zi 

+ He ^ 

+ He^ 


E(MeV/amu) 

Zi 

a(Ex) 

cr(CAL) 

aCAL)/a(Ex) 

2.5 

8 

9.07-20 

7.49-20 

0.826 

2.5 

7 

3.31-20 

3.44-20 

1.04 

2.5 

6 

1.33-20 

1.04-20 

0.782 

2:5 

5 

1.05-20 

5.05-21 

0.481 

1.1 

8 

2.25-18 

2.49-18 

1.11 

1.1 

7 

1.56-18 

1.75-18 

1.12 

1.1 

6 

7.38-19 

7.26-19 

0.984 

1.1 

5 

5.37-19 

4.80-19 

0.894 

1.1 

4 

3.20-19 

2.69-19 

0.840 


*The notation a-b means a x 10'*’. 
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Table 2 

Number of X-Rays per Cosmic Ray Nucleus 


Transition 

X-Ray 

Present Calc. 

Watson (1975) 

Pravdo-Boldt (1975) 

Bussard, et al. (1977) 

Energy (KeV) 

H(a) 

1 (b) 

H 

H 

I 

I 

0+7(2p->Is) 

0.65 

21 

49 (d) 


21.91 

42 . 75 (d) 


0 -*-«( 2 p-ls) 

-0.57 

31 

49 (d) 


18.34 

47 . 21 (e) 


0+6_o+7 

-0.57-0.65 

52 

98 

77(<i) 

40.20 

89.96 


Fe‘*'25 ( 2 p-^ls) 

6.97 

1.4 

e.sO^J) 





Fe’'' 2 '^( 2 p-’-ls) 

~6.5-6.7 

2.3 

9.8(°> 





Pe’’'24_pe+25 

~6.5-6.97 

3.7 

16.3 

11 (d) 



10 . 6 (d) 

Fe+23_Fe+22(3p_2s) 

- 1.1 

43 

117 (d) 





Fe+2i_Fe+i6(3s,3d^2p) 

-0.7-1. 0 

182 

464 





Fe+l6_Fe+23 

-0.7-1. 1 

225 

581 

310(d) 





Calculated for atomic hydrogen medium. 

O’) Calculated for the interstellar medium. 

O’) Includes contributions from charge exchange, excitation and ionization. 
0*^) Contains contributions from excitation and charge exchange. 

(®) Includes charge exchange only. 
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APPENDIX 


Since K- shell fluorescence yields for the ions Fe'*'^ are not available, 

we have fitted the results of theoretical calculations for the "ground state" con- 
figurations of the form Is (2s)^(2p)" . For a given value of n, we have done a 
least squares fit to the function 


OJ. 


(Zp)/ 


1 - ^n(Zp) 


-CZ: 


3 _ 


An + BnZp. 


(Al) 


The constant c was taken to be equal to the value obtained by fitting the data for 
fluorescence yields of neutral atoms with one K-shell vacancy with the form(Al) 
(see Bambynek et al.)» c= - (0,64 ± 0.07) x 10"^ . In(Al) (Zp) is the K-shell 
fluorescence yield for an ion or atom of nuclear charge Zp in the ls(2s)^ (2p)’^ 
configuration. As input data, we have used ajj.’s for neutral B, C, N, and O 
(see Bambynek et al. for references), as well as for various ions: N'*'^ (Z = 1 - 
2) (Bhalla 1975a), Ne*^^ (2 =0^ +5) (Bhalla 1975b, Chen and Crasemann 1975), 
Afi-^ (Z = 3 - 8) (McGuire 1969), and Ar'"^ (Z = 8 - 13) (Bhalla 1973). We note 
that Z indicates the ion’s charge before creation of the K-vacancy. 

The term "ground state" used above means that we assume the projectile 
ion has enough time between collisions to go into the ground state, and further 
more we assume that creation of a K-shell vacancy does not excite the outer 
electrons. For ionization by atomic hydrogen in the interstellar medium, these 
assumptions should be valid. If the target medium becomes denser, and heavier 
species of target atoms become more abundant, the above assumptions break 
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down. In this case, a projectile with a given charge state may contain several 
configurations (Chen and Crasemann 1975), For example Ne"^"^ may contain the 
configurations (Is) ^(2 s)^2p0 , (ls)^2s)^ (2p)^ , and (ls)^(2s)°(2p)2 . The fluores- 
cence yields averaged over the multiplets of each configuration are 0, 0.237, 
and 0.238. After averaging these over the three configurations, assuming sta- 
tistical weighting proportional to the nuniber of multiplets for each configuration, 
a final value of 0.229 is obtained (Chen and Crasemann 1975). Thus in this 
case the ''ground state” result would be coj^ = 0, but the statistically averaged 
result is quite large. 

Using the statistically averaged fluorescence yields for Ne-ions (Chen and 
Crasemann 1975) and for N-ions (Bhalla 1975a), we have again fit the function 
(Al) and determined A^^ and . Since for each n only two data points are avail- 
able, this procedure is not very reliable. 

Table 3 lists the values of the parameters A^^ and of eq. (Al) for the 
charge states arising from the ground-state configurations (Is)^ (2s) ^(2p)" and 
from statistically averaging over possible L-shell configurations. Also given 
are the estimated values of fluorescence yields for Fe-ions obtained by using 
the given A^^ and B^ and letting Zp ^ 26 in (Al). 

In Fig. 16 the ct>j^(26) of Table 3 are plotted as a function of ion charge (be- 
fore the K-vacancy is produced). 

To estimate the effect of the different assumptions concerning the state of 
the ion containing the K-shell vacancy, in Fig. 17 we have plotted the X-ray 
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production cross sections for ~ 6.5 - 6.8KeVX-rays by Fe-ions. In Fig. 17a, 
we give the results Fe-ions in atomic hydrogen assuming ions in the ground 
state just before and after K-vacancy production (dashed line) and assuming the 
statistical distribution of L-shell vacancies after production of the K-shell va- 
cancy (solid line). The multiplicities for these two cases are 2.3 and 2.5, re- 
spectively. In Fig. 17b similar results are obtained for Fe in the interstellar 
medium. In this case the multiplicities are 9.8 (ground state configuration) and 
11.3 (statistical distribution). In all cases the cross sections are calculated as 
in eq. (34) with expression (35) substituted for the last term. 

From Fig. 17 it is seen that the inner-shell Ionization contribution to X-ray 
production can be appreciable in the case of ions with a statistical distribution 
of L-shell vacancies. The correct fluorescence yields to use will depend on the 
assumptions made concerning the target density. 
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Table 3 

The coefficients and Bjj of eq. (Al) and the fluorescence yields obtained for Fe-ions 
when Zp->26 in (Al). The subscript n indicates the number of 2p electrons in the 
“ground state” configuration (IS)^ (2S)^ (2p)". 


Ground' State - Statistically Averaged ■ 

Configuration Configurations 


n 

An 

Bn 

Wn(26) 

An 

Bn 

Wn(26) 

0 

0 

0 

0 

-1.34-1 

8.72-2 

0.953 

1 

2.27-2 

3.61-2 

0.442 

5.87-2 

4.96-2 

0.762 

2 

3.41-2 

3.64-2 

0.468 

2.42-2 

4.25-2 

0.610 

3 

3.75-2 

3.56-2 

0.450 

-3.0 -3 

4.03-2 

0.532 

4 

1.73-2 

3.61-2 

0.443 




5 

3.09-2 

3.50-2 

0.436 




6 

■ 7.46-2 

3.17-2 

0.381 





S6 





FIGURE CAPTIONS 


Figure 1. 


Figure 2. 


Figure 3. 


Figure 4. 


Figure 5. 


Total ionization cross sections for Fe-tons of charge i = +20-+25 
being ionized by H-atoms. The solid curves are the present FBA 
results in cm^ as a function of energy in MeV /nucleon. The dashed 
curves are the BEA results for ionization by protons, and the data 
are from Berkner, et al. (1977), 

Total ionization cross sections for 0~lons of charge i = +3— i-7 being 
ionized by H-atoms. The solid curves are from the FBA calculation 
and the dashed curves are the result of using the BEA for protons as 
targets. 

Total excitation cross sections for the Is 2p transition in 
colliding with H, and He and for colliding with H, Cascading 
from higher excited levels is included as in Equation (22). 

Electron capture by protons from different shells of neutral iron. 

The incident proton is in the energy range 5-40 MeV. It is seen that 
in the few MeV energy range, capture»from the K-shell has the low- 
est cross section. At high energies capture from the K-shell domin- 
ates, as is expected theoretically. 

Cross sections for electron capture by a beam of Fe"^^® from H, He, 


C, N, O, Ne, Mg, Si, S, and Fe, and the total cross section. Cross 
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Figure 6. 


Figure 7. 


Figure 8. 


Figure 9. 


Figure 10. 


section for each, atom A is multiplied by N^/Nh* which is the rela- 
tive interstellar abundance of atom A relative to atomic hydrogen. 
Values of and Njj are taken from Allen (1973). 

Fractions of a beam of ojg^^gen ions in hydrogen and in a mixture of 
gases representing the interstellar medium. The beam fractions 
in charge states to are plotted as a function of energy: 

present calculation, Pravdo-Boldt (1975). The curves 

are labeled by the number of electrons on the ion in that charge 
state. 

Fractions of a beam of Fe-ions in hydrogen and in the interstellar 
medium with charge states Fe"‘'26 to Fe'^^®, as a function of energy. 
The curves are labeled by the number of electrons on the ion in that 
charge state. 

Effective cross section for production of 0.65KeV X-rays by O'^’^ 
in the interstellar medium as a function of ion energy in MeV/amu, 

Effective cross section for production of ~0. 57KeV X-rays by 
O'^^-ions the interstellar mediinn. 

Effective X-ray production cross sections for 0-ions in hydrogen 
as a function of ion energy. 
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Figure 11. 


Figure 12. 


Figure 13. 


Figure 14. 


Figure 15. 


Figure 16. 


Effective cross section for production of 6.97KeV X-rays by Fe"^^^ 
in hydrogen (H) and in the interstellar medium (I), as a function 
of ion energy in MeV/amu, 

Effective cross section for production of ~6,5-6,68KeV X-rays 
by Fe-ions in the interstellar medium. 

Effective cross section for production of ~6.5-6.68KeV X-rays 
by Fe-ions in hydrogen. 

X-ray production cross sections for in (a) the inter- 

stellar medium and (b) hydrogen. 

Effective cross section for production of ~6, 5-6, 8 KeV X-rays by 
Fe-ions in the interstellar medium. This case differs from Figure 
12 in that all capture cross sections are reduced by a factor of 0. 2, 
The multiplicity becomes 5.2, which compares with 9. 8 from 
Figure 12. 

Estimated K-shell fluorescence yields for Fe-ions with charge 
states +16 - +22. The upper curve results from assuming a sta- 
tistical distribution of vacancy states in the L-shell, for a given 
charge state. The lower curve is obtained if the ion remains in the 


"ground” state after a K-shell vacancy is created. 



Fi^re 17. A comparison of the effects of using the two fluorescence yield 
curves of Figure 16 to determine the inner-shell ionization con- 
tribution to the production of ~6. 5 - 6.8KeV X-rays. Part a) is 
the cross section for a hydrogen medium and part b) gives the re 
suits for the interstellar medium. The solid curves correspond 
to the use of fluorescence yields for a statistical distribution of 
vacancy states, while the dashed curves are for the yields from 


a purely "ground" state configuration. 
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